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ABSTRACT 


This  paper  considers  the  possibility  of  inducing  a  convective 

secondary  flow  in  the  fully  developed  channel  flow  of  a  quasi- 

incompressible  (Boussinesq)  fluid.  Instabilities  of  this  type  can 

only  occur  when  the  temperature  gradient  in  the  direction  of  the  body 

force  exceeds  a  certain  critical  value.  This  temperature  gradient  is 

proportional  to  the  Rayleigh  number  of  the  fluid.  We  find  that  for 

channels  of  arbitrary  cross  section,  the  critical  Rayleigh  number  is 

R  ^  1360  (h/d)^  where  h  is  the  arbitrary  channel's  maximum  dimension 
c 

in  the  body  force  direction  and  d  is  the  diameter  of  an  equal  area 
circular  channel.  For  two  special  geometries  it  is  possible  to  im¬ 
prove  the  above  lower  bound  estimate  to  the  critical  Rayleigh  number. 

In  a  circular  channel  R  2:  3450  and  in  a  square  channel  R  ^  2480. 
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INTRODUCTION 

Maslen  [l]  considered  the  fully  developed  steady  flow  of  an  in¬ 
compressible  fluid  in  a  channel  whose  generators  lie  parallel  to  an 
axis.  He  demonstrated  that  the  velocity  components  normal  to  the 
channel  axis  must  vanish  everywhere.  Velte  [2]  treated  the  fully  de¬ 
veloped  steady  flow  of  a  quasi-incompressible  (Boussinesq)  fluid  in  a 
similar  channel  with  a  horizontal  axis.  A  body  force  acts  along  a 
normal  to  the  axis.  The  channel  wall  is  nonuniformly  heated  to  estab¬ 
lish  a  constant  temperature  gradient  in  the  fluid  in  the  direction  of 
the  body  force  ("heated  from  below").  Velte  found  that  there  is  a 
critical  value  for  the  Rayleigh  number  below  which  the  transverse 
velocity  components  vanish.  The  tayleigh  number  is  proportional  to 
the  imposed  temperature  gradient  and  is  given  by  R  ■  gaf3h  /kv  where  a 
is  the  fluid  coefficient  of  thermal  expansion,  g  is  the  body  force  per 
unit  mass  acting  on  the  fluid,  j3  is  the  magnitude  of  the  imposed  tem¬ 
perature  gradient,  h  is  the  maximum  dimension  of  the  region  in  the 
direction  of  the  body  force,  k  is  the  fluid  thermal  diffusivity,  and 
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v  is  the  fluid  kinematic  viscosity.  Above  the  critical  Rayleigh  num¬ 
ber  it  is  possible  to  establish  a  convective  secondary  flow. 

Velte  developed  a  method  for  calculating  upper  bounds  to  this 
critical  Rayleigh  number  in  channels  of  arbitrary  cross  section.  If 
for  design  purposes  it  is  desirable  to  Insure  the  stability  of  the 
initial  fully  developed  channel  flow,  then  lower  bounds  to  the  criti¬ 
cal  Rayleigh  number  must  be  established.  In  this  paper  a  method  is 
presented  to  calculate  these  lower  bounds. 
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THE  GOVERNING  EQUATIONS 


Consider  a  fully  developed  quasi-incompressible  flow  in  a  channel 
parallel  to  the  y-axis.  The  Cartesian  velocity  components  u,  v,  and  w 
are  in  the  x,  y,  and  z  directions  (given  by  the  unit  vectors  i^,  and 
k).  A  body  force  g  acts  in  the  negative  z  direction.  The  initial 
velocity,  temperature,  and  pressure  distributions  in  the  fluid  are 


u  «  0,  v  -  v  (x,  z)j,  w  «  0 

o  l 


VT  ®  -0k 


Vp 


-pg(l  -  af3z)k  +  pvv  v  j 


> 
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where  T  is  the  temperature,  p  is  the  pressure,  p  is  the  mean  fluid 

2 

density,  V  is  the  gradient  operator,  and  V  is  the  Laplacian  operator. 

It  is  assumed  that  the  secondary  convective  flow  is  fully  develop¬ 
ed  (independent  of  y),  hence  the  governing  dimensionless  perturbation 
equations  valid  at  the  onset  of  instability  are  (see  [2]) 


vViIr  =  R0 


(2) 
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(3) 


where  il  is  the  perturbation  stream  function  measured  in  units  of  k 

(u*-t,w»i|r)#  and  0  is  the  temperature  perturbation  in  units  of 
%  X 

Sh.  Here  the  gradient  and  Laplacian  operators  are  ir.  the  x,  z-plane. 
If  the  rigid  channel  wall  has  a  large  thermal  conductivity  and  heat 
capacity  relative  to  the  fluid,  then  the  initial  temperature  distri¬ 
bution  on  the  wall  is  maintained  at  all  times.  Thus  the  boundary  con- 
di.  i.ons  on  the  channel  wall  B  are 


i!f  ■  Vilf  *  n  *  0  on  B 

0*0  on  B 


(4) 


Here  n  is  a  unit  outer  normal  to  B.  The  above  eigenvalue  problem  is 
equivalent  to  the  variational  principle  in  the  cross-section  plane  C. 

/  (V2t)2  dx  dz 

R  *  ■  .c - -  ■■  — - ,  i)f*Vf'n*0*O  on  B  (5) 

JC(V0)2  dx  dz 

The  functions  i|f  and  0  are  chosen  from  a  space  of  admict4ble  func¬ 
tions  whose  f  functions  are  four  times  differentiable  and  whose  0 

functions  are  twice  differentiable.  In  addition,  the  functions  satis- 

2 

fy  the  constraint  V  0  ■  i|f  in  the  region  C  and  the  boundary  conditions 

X 

on  B.  The  critical  Rayleigh  number  is  given  by  R  *  min  R. 

c 

GENERAL  LOWER  BOUND  ESTIMATES  TO  R 

L  1  —  —  ~  . .  ■— c 

2  2  2 

Using  the  inequality  (V<f)  ^  f  »  and  the  constraint  V  0  ■  t  ,  it 

X  X 

follows  directly  from  equation  (5)  that 
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J  (V2*)2  dx  dss  J  (7Z0)Z  dx  dz 

R  ^  min  — — - -r -  *  min  — - 5 - 

c  J*  (V\|f)  dx  dz  J  (V6)  dx  dz 


2rtv2 


over  the  function  space  2^.  Employing  the  Principle  of  Monotony  [3], 

one  finds  that  a  lower  bound  to  R  is  also  given  by  the  extremum  prob- 

c 

lem  (6)  over  a  less  restrictive  function  space  2^  *n  which  f  and  9  are 
chosen  independent  of  one  another.  Hence  2^  is  a  sub-space  of  22* 
Consider  the  extremum  problem 


Jc(vV  dx  dz 
J  (V^)2  dx  dz 


i|r  «  V\|f  •  n  *  0  on  B  (7) 


over  the  function  space  2^.  The  smallest  extremal  value  of  equation 
(7)  is  identical  to  the  principal  eigenvalue  of  the  problem  posed  by 

V2?2*  +  sVijr  -  0  in  C 


t  ■  Vtjf  •  n  =*  0  on  B 


For  most  regions,  the  principal  eigenvalue  cannot  be  calculated  in  a 
simple  manner.  But  for  a  circular  channel  (in  which  h  «  d,  the  diame¬ 
ter),  K2  is  determined  from  the  first  root  of  J^(K/2)  *  0  (see  [4]) 
where  is  the  Bessel  function  of  the  first  kind  of  order  one.  There- 
fore,  k  «  58.74.  Polya  and  Szego  [4]  have  demonstrated  that  among 

all  plane  regions  of  equal  area,  the  circle  has  the  smallest  value  of 
2 

K  .  Thus  if  one  replaces  a  given  arbitrary  region  of  height  h  by  an 
equal  area  circle  of  diameter  d,  then 


K2  *  58.74  (h/d)2 


Next  consider  the  extremum  problem 


2  J*  (V9)2  dx  dz 

IT  »  min  — —  -r - - - ,  0  *  0  on  B  (10) 

f  6  dx  dz 
d  c 

over  the  function  space  The  smallest  extremal  vaiue  of  equation 
(10)  is  identical  to  the  principal  eigenvalue  of  the  following  prob¬ 
lem: 


V20  +  L20  -  0  in  C 

(ID 

0  ■  0  on  B 

2 

For  a  circle,  L  is  .^termined  from  the  first  root  of  J  (L/2)  ®  0. 

2  * 

Therefore,  h  -  23.14.  Polya  and  Szego  also  have  shown  that  among 

all  plane  regions  of  equal  area,  the  circle  has  the  smallest  value  of 
2 

L  .  Analogous  to  the  estimate  (9),  one  then  can  obtain 


L2  *  23.14  (h/d)2 


(12) 


Using  the  divergence  theorem,  the  boundary  condition  8  =  0  on  B,  and 
the  Schwarz  inequality,  it  follows  that 


ir 


0V20 


dx  dz 


2 


£ 


dx  dz 


dx  dz 


Combining  these  expressions  with  equation  (10)  yields 
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f  (v2e>2  dx  dz 
Jc  (V0)2  dx  dz 


(13) 


Substituting  the  estimates  (7)  and  (13)  into  equation  (6)  and  using 
the  estimates  (9)  and  (12)  gives 


Rc  *  1360  (h/d)4 


(14) 


As  an  illustrative  example,  consider  a  channel  whose  cross  section 
is  an  equilateral  triangle  of  height  h.  A  simple  calculation  yields 

A 

(h/d)  *  1.85  and  R£  2s  2520  for  this  channel. 

LOWER  BOUNDS  TO  P  IN  TWO  SPECIAL  CASES 

For  channels  with  circular  or  square  cross  sections,  the  lower 
bound  estimates  of  the  previous  section  may  be  improved.  In  these  two 
special  cases  we  find  estimates  from  equation  (6)  over  the  function 
space  rather  than  the  less  restrictive  function  space 

The  Circular  Channel.  In  a  circular  channel  (0  s:  r  £  %, 

0  £  cp  £  2n,  cp  measured  counterclockwise  from  the  x-axis)  the  complete 
set  of  eigenfunctions  to  equation  (8)  is  given  by  (see  ref.  [4]) 


eo  oo 


l  I  {<2r>m  JJcf) "  4(Kra„0} 


to«*0  n»l 


iA  cos  top  +  B  sin  mcpr 
1  ran  ^  ran  v 


(15) 


2  fch 

The  associated  eigenvalues  are  K  where  K  is  the  n  root  of 

um  ran 

Jatfl(K/2)  *  ^el-e  [2]  has  demonstrated  that  the  stream  function 
mode  associated  with  the  critical  Rayleigh  number  is  synmetric  with 


respect  to  both  the  x-  and  the  z-axis.  Thus  the  subset  of  (15)  that 
has  the  proper  syramecry  property  is  given  by 


*  •  l  l  AJ<2r>m  Jm(rf)  •  Jm  fan1)}  c°s 

m=>0  n*l 

m  an  even  integer  (16) 

2  2 

The  principal  eigenvalue  of  the  subset  is  K  «  *  58.74* 

The  complete  set  of  eigenfunctions  to  equation  (11)  for  the  cir¬ 
cular  channel  is  given  by 

00  00 

6  -  1  I  7mOmnO  {Amn  cos  ®mn  sin  (17) 

m=*0  n«l 

2  th 

The  associated  eigenvalues  are  L  where  L  is  the  n  root  of 

mn  mn 

2 

Jm(L/2)  «  0.  From  the  constraint  V  0  «  i|f  and  the  synmetry  property 
of  ilr ,  one  can  conclude  that  0  must  be  symmetric  with  respect  to  the 
x-axis  and  anti -symmetric  with  respect  to  the  z-axis.  The  complete 
subset  of  (17)  that  has  the  proper  symmetry  property  is  given  by 


00  00 

9  ”  1  Z  AmnJmCLmnO  cos  mcp,  m  an  odd  integer  (18) 
m»l  n»l 

2  2 

The  principal  eigenvalue  of  the  subset  is  L  ■  ■  58.74.  One  now 

has  at  hand  the  following  estimates  for  the  circular  channel 
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0  L  (V2'!')2  dx  dz 

K*  *  — - x - 

01  p 


J  (V<0  dx  dz 


9  Sr  OA)2  dX  dZ 

T  “  <  .9. .  . .  ■ 

L11  r  /„-v2 


f  (Vi)  dx  dz 

J  c 


Substituting  the  above  estimates  into  equation  (6)  gives 


R  *  3450 

c 


(19) 


for  the  circular  channel.  Recently  Sherman  [5]  used  a  Rayleigh-Ritz 
technique  and  the  variational  principle  (5)  to  establish  the  upper 
bound  estimate 


Rc  £  6510  (20) 

for  the  circular  channel. 

The  .Square  Channel-  In  a  square  channel  (jx|  £  %,  |z|  *  h),  the 

eigenfunctions  of  equation  (3)  cannot  be  determined  in  any  simple 

form.  Weinstein  [6]  has  determined  a  lower  bound  to  the  principal 
2  2 

eigenvalue,  K  *:  5.1  tt  .  The  associated  principal  eigenfunction  is 
symmetric  with  respect  to  the  x-  and  the  z-axis.  Velte  [2]  has  shown 
that  the  critical  stream  function  mode  also  has  this  symmetry. 

The  complete  set  of  eigenfunctions  tc  equation  (11)  for  the 


square  channel  is 
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1  1  K cos  c2n 


1)  ttx  +  B  sin  2mTTxr 
m  j 


m»l  n«l 


^An  cos  (2n  -  1)  rry  +  Bn  sin  2«TT/j- 


As  in  the  circular  channel,  one  can  conclude  that  9  must  be  symmetric 
with  respect  to  the  x-axis  and  anti -synme trie  with  respect  to  the  z- 
axis.  The  complete  subset  of  eigenfunctions  with  the  required  symme¬ 
try  is 


l  l 

m-1  n-1 


A  sin  2mTTx  cos  (2n  -  1)  ny 
mn  \  /  j 


2  2 

The  principal  eigenvalue  of  this  subset  is  L  ■  5tt  .  Thus  for  the 
square  channel  one  has  the  estimates 


,  L  <*2t>2  <ix  i* 

5.1  «  •? - 5 - 

L  Mr  dz 


Jc  (vV  dx  dz 

^  n 

jc  (?0)  dx  dz 


Substituting  the  above  estimates  into  equation  (6)  gives 


R  *  2480 
c 


Velte  calculated  the  upper  bound  estimate  for  the  square  channel 


R  *  5030 
c 


(24) 


CONCLUDING  REMARKS 


The  lower  bounds  to  established  in  this  paper  differ  con¬ 
siderably  from  the  upper  bounds  calculated  by  the  Rayleigh-Ritz 
method.  Nevertheless  the  lower  bounds  provide  the  designer  with  a 
limit  of  error  which  the  Rayleigh-Ritz  method  cannot  yield.  Better 
lower  bounds  may  be  established  if  the  Weinstein  method  [3]  can  be 
successfully  applied  to  the  variational  principle  given  by  equation 
(5).  It  is  evident  from  both  lower  bound  and  upper  bound  calculations 
that  the  nature  of  the  confining  region  can  have  a  marked  effect  on 
the  critical  Rayleigh  number. 

Finally  it  should  be  noted  that  for  rully  developed  flow  between 
two  horizontal  planes  of  infinite  lateral  extent,  the  critical  Rayleigh 
number  is  exactly  R£  *  1707.8  [7]» 
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